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Abstract 

The phase transition temperature for the Bose-Einstein condensation (BEC) 
of weakly-interacting Bose gases in three dimensions is known to be related to 
certain non-universal properties of the phase transition of three-dimensional 
0(2) symmetric c/)^ theory. These properties have been measured previously 
in Monte Carlo lattice simulations. They have also been approximated analyt- 
ically, with moderate success, by large N approximations to 0{N) symmetric 
(j)^ theory. To begin investigating the region of validity of the large N approx- 
imation in this application, I have applied the same Monte Carlo technique 
developed for the 0(2) model (Ref. [5]) to 0(1) and 0(4) theories. My results 
indicate that there might exist some theoretically unanticipated systematic 
errors in the extrapolation of the continuum value from lattice Monte Carlo re- 
sults. The final results show that the difference between simulations and NLO 
large A'^ calculations does not improve significantly from N = 2 to N = 4:. 
This suggests one would need to simulate yet larger A^'s to see true large 
scaling of the difference. Quite unexpectedly (and presumably accidentally), 
my Monte Carlo result for = 1 seems to give the best agreement with the 
large A'^ approximation among the three cases. 

I. INTRODUCTION 

The computation of the phase transition temperature Tc for dilute or weakly-interacting 
bose gases has attracted considerable interest. Due to the non-perturbative nature of long- 
distance fluctuations at the second-order phase transition at large distance, the calculation 
of corrections to the ideal gas formula for Tc is non-trivial. In the dilute or weakly-interacting 
limit, the correction AT^ = — Tq to the ideal gas result Tq for a homogeneous gas can be 
parametrized as^ 

AT 

^ = ca,,nV3 + 0(aln'/' Ha^^n'/')), (1.1) 
-t n 



^For a clean argument that the first correction is linear in Osc, see Ref. [1]. For a discussion of higher-order 
corrections, see Refs. [2,3]. 



1 



where age is the scattering length of the 2-particlc interaction, n is the number density of the 
homogeneous gas. c is a numerical constant, and 0(- ■ ■) shows the parametric size of higher- 
order corrections. Baym et al. [1] have shown that the computation of c can be reduced 
to a problem in three dimensional 0(2) 0^ field theory. In general, O(A^) (j)^ field theory is 
described by the continuum action 



(1.2) 



where = (0i, 02) ■ ■ ■ j 4'n) is an A^-component real field. 1 will focus exclusively on the case 
where r has been adjusted to be at the order/disorder phase transition for this theory for a 
given value of the quartic couphng u. The relationship to for Bose-Einstein condensation 
found by Baym et al. is that the constant c in (1.1) is^ 
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where 0^ = 0f + 0^ H (f)% , N ^ 2 and 



A(0^ 



(0^ 
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represents the difference between the critical-point value of (0^) for the cases of (i) u non- 
zero and (ii) the ideal gas u—0. Thus the computation of the first correction to Tc due to 
interactions is related to the evaluation of A(0^)c in three-dimensional 0^ theory. 

Having tuned r to the phase transition, u is then the single remaining parameter of the 
three-dimensional continuum 0^ theory (1.2). The dependence of (ultraviolet-convergent) 
quantities on u is determined by simple dimensional analysis, and u has dimensions of inverse 
length. The A(0^)c/it in (1.3) is dimensionless and so is a number independent of u in the 
continuum theory. Monte Carlo simulations of this quantity in 0(2) theory have given 
c = 1.29 ± 0.05 [4] and c = 1.32 ± 0.02 [5]. 

One of the few moderately successful attempts to approximate this result with an analytic 
calculation has been through the use of the large N approximation. (But see also the 4th- 
order hnear S expansion results of Refs. [7] . For a brief comparison of the results of a wide 
spread of attempts to estimate c, see the introduction to Ref. [6] and also Ref. [8].). The 
procedure is to calculate A(0^)c/'U for O(A^) theory in the limit that N is arbitrarily large, 
and then substitute the actual value N=2 of interest into the result. The large result 
was originally computed at leading order in 1/N by Baym, Blaizot, and Zinn-Justin [9] and 
extended to next-to-leading order (NLO) in 1/N by Arnold and Tomasik [10], giving 



A(02 



N 



u 



967r2 



(1.5) 



^This is given separately in Ref. [1] as ATc/Tq = —2mkBTo A{(jP)c/3fi^n and the identification of u as 
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-0.000494(41) 


-0.0004990 


- 1(8) % 


0.0015249(48) 


2 (Ref. [5]) 


-0.001198(17) 


-0.001551 


+30(2) <X 


0.0019201(21) 


4 


-0.00289(18) 


-0.003665 


+27(8) % 


0.002558(16) 



TABLE I. Simulation results for 7V=1,2,4 and the NLO large N results for A{(j)'^)c/u. The difference 
column shows the percentage excess of the magnitude of the large N approximation result over the magnitude 
of the simulation result. N = 2 simulation results are quoted from Ref. [5]. 

Setting N — 2, one obtains c ~ 1.71, which is roughly 30% higher than the results obtained 
by Monte Carlo simulation of 0(2) theory. Considering that N — 2 has been treated as 
large in this approximation, the result is fairly encouraging. 

The goal of the present work is to further explore the applicability of the large result 
(1.5) by testing it for other values of N. I have applied to other 0(iV) models the same 
techniques used in Ref. [5] to simulate the 0(2) model. In this paper I report the measure- 
ment of A(0^)c/m for the 0(1) and 0(4) theories. The final results, compared to the large 
N approximations of (1.5), arc given in table I. 

As a byproduct of the analysis, 1 also report the measurement of the critical value of 
r. The coefficient r requires ultraviolet renormalization and so is convention dependent. In 
Table I, I report the dimensionless continuum values of Vc/u'^ with Vc defined by dimensional 
regularization and modified minimal subtraction (MS) renormalization at a renormalization 
scale /i set to m/3. Among other things, this quantity can be related to the coefficient of the 
second-order (a^^n^^^) correction to ATc [2,11]. One can convert to other choices of the MS 
renormalization scale p, by the (exact) identity 

V? 18(47r)2 //2 

In section II, I briefiy review the algorithm and the necessary background materials and 
formulas for improving the extrapolations of the continuum and infinite- volume limits. Most 
of the technical details can be found in Ref. [5]. Section III gives the simulation results and 
analysis. Section IV is the conclusion from comparing the numerical results with the NLO 
large N calculation (1.5). 

II. LATTICE ACTION AND METHODS 

The action of the theory on the lattice is given by 

Slat = a' E{^*iat(-VL)^iat + \roK, + ^(^D^, (2-1) 

where a is the lattice size (not to be confused with the scattering length a^^)- I will work on 
simple cubic lattices with cubic total volumes and periodic boundary conditions. As in Ref. 
[5] and as reviewed further below, the bare lattice operators and couplings ($iat, ''^o, 
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(1.6) 



Uo) are matched to the continuum, using results from lattice perturbation theory to improve 
the approach to the continuum limit for finite but small lattice spacing. 
I shall use the improved lattice Laplacian 

V'$(x) = a-' Yl [-T2^i^ + + + ai) - |$(x) + |$(x - ai) - ^^^(x - 2ai)] , 

i 

(2.2) 

which (by itself) has O(a^) errors, rather than the standard unimproved Laplacian 

Vu$(x) = Yl [*(^ + «i) - 2$(x) + *(x - ai)] , (2.3) 

i 

which has O(a^) errors. Readers interested in simulation results with the unimproved kinetic 
term should see Ref. [23]. 

Because the only parameter of the continuum theory (1.2) at its phase transition is u, the 
relevant distance scale for the physics of interest is of order 1/uhy dimensional analysis. To 
approximate the continuum infinite- volume limit, one needs lattices with total linear size L 
large compared to this scale and lattice spacing a small compared to this scale. As a result, 
ua is the relevant dimensionless expansion parameter for perturbative matching calculations 
intended to eliminate lattice spacing errors to some order in a. The analogous dimensionless 
parameter that will appear in the later discussions of the large volume limit is Lu. 



A. Matching to continuum parameters 

To account for the lattice spacing errors, I have adapted the perturbative matching 
calculations to improve lattice spacing errors given in Ref. [5], where the details of the 
calculations can be found. Here I will simply collect the results from that reference. The 
lattice action expressed in terms of continuum parameters is written in the form 



'lat 



El^(Viat</')^ + §^(r + <5r)0^ + 



(2.4) 



The continuum approximate value for A(0^)c is theoretically expected to have the form 



A((/.2)e = Z.(</.')lat-<^</'' + 0(a2). 



(2.5) 



For the improved Laplacian, the matching calculations have been done to two-loop level (2/) 
for Z^, Zr, Sr and 6u, and three loops (3/) for 50^, yielding 
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(2.6b) 
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6 An"" ■ 18(47r)2 
where, for the improved Laplacian (2.2 



Ca + ln 
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(47r)2 



ma 



E ~ 2.75238391130752, 
{ ~ -0.083647053040968, 

Ci - 0.0550612, 

C2 - 0.03344161, 

C3 ~ -0.86147916, 

C4 ~ 0.282. 



(2.6c) 
(2.6d) 
(2.6e) 



(2.7a) 
(2.7b) 
(2.7c) 
(2.7d) 
(2.7e) 
(2.7f) 



As shown in Ref . [5] , the result of this improvement is that at fixed physical system size 
Lu, the lattice spacing error of A{(f)'^)c/u should be 0(0^). However as will be shown in the 
following, my simulation results indicate that there might still exist some linear coefficients 
even after applying the formula (2.5). I assign an error to my continuum extrapolations that 
covers both linear and quadratic extrapolations. 

This has produced about 10% systematic error in the final value of A{(f)'^)c/u. The 
extrapolation of Tc/m^ on the other hand has 0{a) error, since no improvement has been 
made. Overall, the fitting formulas for the data taken at fixed Lu are given as 



u 



Lu 



U^JLu 



Bi + B2{uaf, 
Di + D2{ua), 



(2.8a) 
(2.8b) 



for a quadratic fit of A(02)c/w, and 



u 



Lu 



Lu 



B, + B!,(ua), 
Di + D2{ua), 



(2.9a) 
(2.9b) 



for a linear fit of A{(j)'^)c/u. 



B. Algorithm and Finite Volume Scaling 

Working in lattice units (a=l) with the lattice action (2.4), I update the system by 
heatbath and multi-grid methods [21]. At each level of the multi-grid update, I perform over- 
relaxation updates. Statistical errors are computed using the standard jackknife method. 
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The strategy is to vary r at fixed u to reach the phase transition point. In order to define 
a nominal "phase transition" point in finite vohimc, I use the method of Binder cumulants 
[22] . The Binder cumulant of interest is defined by 



where is the volume average of 0, 

0=i| d'x(t>{x). (2.11) 

The nominal phase transition is defined as occurring when C = C*, where C* is a universal 
value which improves convergence to the infinite volume limit. C* depends on the shape 
and boundary conditions of the total lattice volume but not on the short distance structure. 
For cubic lattice volumes with periodic boundary conditions, C* ~ 1.603(1) for 0(1) [16] 
and C* ~ 1.095(1) for 0(4) [17]. I have checked that the errors on C* are not significant for 
the purposes of this application. Therefore I have used the central values for the nominal 
C*. 

In practice, it never happens that the simulation is done precisely at the r where C — C* . 
I instead simply simulate at some r = rgim close to it. I then use canonical reweighting of 

the time series data to analyze r's near rgim and determine r and A(0^) at C = C* . 

A renormalization group analysis of the scaling of finite volume errors shows that, when 
the method of Binder cumulants is used, the values of A(0^)c and scale at large volumes 
as [5] 

~ A + A^iLuy+y^ (2.12a) 

u 



^ c^Ci + C2{Lu)-y'-^ (2.12b) 



for fixed ua. Here d—3 is the dimension of space, yt = l/i', and u and u are the standard 
0(A'") critical exponents associated with the correlation length and corrections to scaling, 
respectively. The values of the exponents that I have used are^ 

0(1): yt = 1.587(1), 

^ = 0.84(4), 

0(4): yt - 1.329(2), ^^•'"'^ 

cu = 0.79(4). 



•^A nice review and summary of the critical exponents can be fomid in [12]. For the case of 0(1), we have 
used their summarized values of = l/jji and lu based on the calculation from High Temperature (HT) 
expansion technique and Monte Carlo simulations. For comparison, some experimental results for yt are: 
1.61(8) [13], 1.58(2) [14]. For 0(4), Monte Carlo simulation for gives: 0.7525(10) [15](used by us), 0.749(2) 
[18]. From HT expansion: 0.759(3) [19]; from e-expansion 0.737(8) [20]. For uj, the only MC simulation 
result is w = 0.765 (without quotation of error) from [18]. Ref. [20] gives u) = 0.774(20) (d=3 exp.) and 
0.795(30) (e-exp.). We have chosen the value to cover both of the errors for u. 
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The large volume scaling of A{(f)'^)c/u depends only on yt, which have errors ~ 0.1% in both 
cases. They have a negligible effect on my final results. On the other hand, for the case of 
Tc/u^, the large volume scaling depends on both yt and cu. I will show that the errors of lu 
do have effects on the large volume errors of rdu^. They are included in the final results. 

For a discussion of higher-order terms in the large volume expansion, see Ref. [5], but 
the above will be adequate for this application. I will fit the large Lu data to the leading 
scaling forms (2.12) to extract the infinite volume limit. 

The basic procedure for extracting simultaneously the wa — > and Lu ^ oo limits of 
my result will be as follows, (i) I first fix a reasonably small value of (ua)* of ua, take data 
for a variety of sizes Lu (to as large Lu as practical) , and then extrapolate the size of finite 
volume corrections, fitting the coefficients A2 and C2 of the scaling law (2.12). (ii) Next, I 
instead fix a reasonably large physical size (Lu)* and take data for a variety of ua (to as 
small ua as practical), and extrapolate the continuum limit of our results at that {Lu)*, 
which corresponds to fitting Bi and Di of (2.8, 2.9). (iii) Finally, 1 apply to the continuum 
result of step (ii) the finite volume correction for (Lu)* determined by step (i). In total, I 
have 



u 



> =B^-A2{{Luy)-''+y\ (2.14) 

final 



J final 



- C2((L«)*)-^*-'^. 



There is a finite lattice spacing error in the extraction of the large volume correction {i.e. 
A2 and C2). In Ref. [5], it is argued that this source of error is formally high order in {ua)* 
and so expected to be small. 



III. SIMULATION RESULTS 
A. 0(1) Results for A{(^^)c/u 

There is a practical tradeoff between how large one can take the system size {Lu)* in 
order to reach the large volume limit versus how small one can take {ua)* in order to reach 
the continuum limit. Fig. 1(a) and the lower points (circles) in 2(a) show, respectively, the 
Lu dependence at {ua)* = 12 and the ua dependence at {Lu)* = 768 for A{(l)'^)c/u in the 
0(1) model. From the Lu dependence, we can see that {Lu)* = 768 is a reasonably large 
value of Lu — the finite volume corrections are moderately small. From the ua dependence, 
we can see that {ua)* = 12 is reasonably small. 

Fig. 1(b) shows the size A2 x {Lu)~'^~y* at {Lu)* = 768 when fitting the Lu dependence of 
Fig. 1(a) to the scaling form (2.12a). The result of the fit, and its associated confidence level, 
depends on how many points are included in the fit. Percentage confidence levels are shown 
in the figure. My procedure will be to determine the best values of the fit parameters by 
including as many points as possible while maintaining a reasonable confidence level. Then 
to assign an error, I use the statistical error from including one less point in the fit. This will 
help avoid the underestimation of systematic errors. The resulting estimate 0.0001335(50) 
of the finite volume correction is depicted by the shaded bar in Fig. 1(b) and collected in 
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Table II. The corresponding best fit is shown as the solid line in Fig. 1(a). Even though 
there is no direct use of it, I show for completeness the fit parameter Ai in Fig. 1(c), which 
corresponds to the infinite-volume value of A(0^)c/m at (ua)* = 12. 

Having found the large volume extrapolated values, I will now move on to the — > 
hmit for A{(f)^)c/u at a fixed physical volume (Lu)* — 768. Theoretically, from the argument 
of Ref. [5], if one uses the 0(a) improved formula (2.5) to obtain the continuum approximated 
A{(j)'^)c/u, the remaining lattice error should be ~ 0{uay. 

The circular data in Fig. 2(a) show the 0(a) improved simulation results by using Eq. 
(2.5) . First I have tried a quadratic fit with the form Bi + B2{uaf . The fitting results for 
El are given by Fig. 2(b). One can see while a three-point {ua < 8) fit gives an impressive 
C.L. of 85%, including the fourth point {ua = 9.6) reduces the CL. to 17%. However, adding 
two more points keeps the C.L. above 10%. Due to this feature, it is hard to determine the 
last point that should be fit to the quadratic formula. Instead, I assigned the value to cover 
all the 85%, 17%, 13% and 10% C.L.s which gives -0.000383(17). The result is indicated by 
the shaded area. 

The poor behavior of the quadratic fit makes one wonder if the data is actually more 
"linear" than "quadratic." To check this, I re- fit the circular data using a linear function 
(2.9a). The results are shown in Fig. 2(c). Surprisingly, I can fit all the data points with 
a very good C.L.. The ua ^ Q extrapolated value in this case is -0.0003275(55), which is 
very different from the quadratic fit result. The obvious hnear behavior of the data seem to 
indicate that there might be some residual 0(a) coefficient in Eq. (2.5) that has not been 
accounted for. Given this uncertainty, I have used both the linear and quadratic extrapolated 
results, combined with the error due to the finite lattice size to obtain two continuum values 
(Table II). They differ by about 10%, which is considered my systematic error. The final 
result is assigned to cover both results and is tabulated in Tables II and I.^ 

As a comparison, I have also shown the naive subtracted result given by A($^)l/mo = 
((<I>^) — NH/Att) /uq vs. the unimproved u^a as the scale based on the same simulation (the 
stars in Fig. 2(a)). One can see that the finite lattice spacing errors is more explicit with 
the stars, corresponding to the unimproved data, than that of the circles, corresponding to 
the improved ones (using (2.5)). Theoretically A{^^)i^/uo data have log x linear and linear 
dependences on uga. 

B. Tc and 0(4) Results 

Figs. 3 and 4 show similar curves for the 0(4) model, where we have taken (Lu)* = 144 
and (ua)* = 3 as a reasonably large size and reasonably small lattice spacing. It's worth 
noting the large N limit predicts that the distance scale that characterizes the physics of 
interest should scale as 1/ (Nu) . That leads one to expect that the upper limit for "reasonably 



^In Ref. [5], for the case of 0(2) model, the Lu fixed data are fitted by quadratic functions only. The 
Lu = 144 data can be fitted also by a linear function. The Lu = 576 data on the other hand can be fitted 

by only a quadratic function for the first four points. However if one takes off the smallest ua point, the rest 
four points can be fitted very well by a straight line too. This might indicate that there could be also some 
linear coefficients in the 0(2) theory. For more discussions on this possible linear coefficient see Ref. [23]. 
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small" values of ua and lower limit for "reasonably large" values of Lu should scale roughly 

as 1/N. 

The fitting proceeds as in the 0(1) case. For the Lu-Hxed data of A(0^)c/m (Fig. 4), I 
have again fitted the data with both linear and quadratic functions. Both the quadratic fit 
and the linear fit can go up to = 8 while keeping good C.L.s. However their ua ^ 
extrapolations are also different (see Fig. 4 and Table II for the fitting results). The final 
continuum value is assigned to cover both extrapolations. 

The analysis of the result for Tc/u^ in both the 0(1) and 0(4) models is much the same 
(Figs. 5-7). For the large volume extrapolations, since the critical exponent cu has a larger 
uncertainty (~ 5%), I have used three different values of cu (maximum, central, minimum) 
for both 0(1) and 0(4)'s cases. The final results of the C2 x (Lm)^^*^"^ value cover all the 
three different extrapolated values. This difference results in about 1% of the error quoted 
in the final results for Tc/u^. For the Lu fixed data, since there is no improvement of 0(a) 
errors, the continuum extrapolation of the data taken at fixed Lu are fitted to linear behavior 
in ua, following (2.8b or 2.9b). The fitting results are again summarized in Tables I and II. 
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FIG. 1. (a) shows the simulation result for A((/)^)c/u vs. physical volume at ua = 12 for 0(1); (b) is 
the results of the difference between the (768, 12)* point and the infinite extrapolation, (c) is the Lu — > oo 
extrapolations; and Since we used the fitting formula for A{(jP)c/u at (768, 12)* 2.5, this difference depends 
only on the results of A2 but not on Ai. In both (b) and (c), the shaded areas are the quoted results. 
The confidence levels are given as as the percentage numbers on top the fitted values. "NA" stands for 
"non-applicable", which means the number of the fitting parameters equals the number of fitted points. 
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FIG. 2. Simulation results for A((/)^)c/u vs ua at fixed physical volume {Lu}* — 768 for 0(1). The 
circular data are obtained from (2.5), with S(j)^ given by (2.6a). The star data are obtained using the naive 
lattice result A($^)l = {{^^) — N11/A'!t)/uq. It has linear+linearx log+quadratic dependence on ua. (b) 
gives the fitting results for the quadratic fit and (c) gives the fitting results for the linear fit. 
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FIG. 3. As Fig. 1 but for 0(4) at ua = 3 with reference point {Lu,ua)* = (144,3). 
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FIG. 4. As Fig. 2, but for 0(4) at {Lu)* — 144. The triangular data in (a) is again the naive lattice 
result A((f>^)L/uo which shows again bigger lattice errors than the 0{a) improved data(circular data). 
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FIG. 6. As Fig. 5 but for 0(4) at ua = 3 with (Lu,Ma)* = (144,3). 
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N 




^1 + ^2 X {Luy^+y' fit 


Bi + B2{uaf fit 


Bi + B'^{ua) fit 


A(02) J-u = - ^2 X (Lm)-'^+2'* 








0(Ma)2 


0{ua) 


final 


1 


(768,12) 


-0.0006199(40) 


0.0001335(50) 


-0.000383(17) 


-0.003275(55) 


-0.000517(18) 


-0.0004610(74) 


-0.000494(41) 


4 


(144,3) 


-0.003023(92) 


0.002978(92) 


0.000007(26) 


0.000163(37) 


-0.002971(96) 


-0.002815(99) 


-0.00289(18) 




N 


(Lm, mo)* 


+ C2 X (Lm)-?/*-^ fit 


Di + D2{ua) fit 


rju^ = Di - C2 X (Lw)-?'*-'^ 


Ci 


C2 X {Ui)-y'-^ 




1 


(768,12) 


0.00274001(28) 


-0.00000178(22) 


0.0015231(48) 


0.0015249(48) 


4 


(144,3) 


0.0035035(62) 


-0.000183(12) 


0.002375(10) 


0.002558(16) 



TABLE II. Fitting results for A(0^)c/w and Vc/v? using formula (2.15). In the case of A(0^)c/w, both a linear fit result and a quadratic fit result 
Bi are given. The final values are assigned to cover both fit results. 



IV. CONCLUSION 



Fig. 8 and Table II show a comparison of the simulation results for A(0^)c/m with the 
NLO result (1.5) for N — 1,2,4. The N — 1,4 cases are from this paper. The result of 
0(2) model is taken from Ref. [5].^ In the figure, I have actually plotted A{(f)'^)c/{Nu), 
where the explicit factor of N factors out the leading-order dependence on as iV — > oo. 
Amusingly, the large estimate for A^=l seems to be the most accurate of the three cases. 
This is presumably an accident — the result of the large approximation just happens to be 
crossing the set of actual values near N — 1. Prom Eq. (1.5), the error of the for A{(f>'^)c/Nu 
should scale as but there is clearly no sign of such behavior between A^=2 and N—A. 

This might indicate that N=4 is still too small for the error in the large approximation 
to scale properly. It would be interesting to numerically study yet higher A^ such as A^=8 
or A^=16 to attempt to verify the details of the approach to the large A^ limit. 

For the systematic error due to the possible hnear behavior in the — > extrapolation 
of A(0^)c/ti, so far no theoretical reasoning is found. It would be interesting to reinvestigate 
the matching calculations given in Ref. [5]. 
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^If one also considers the systematic uncertainty in the ua — > extrapolation in Ref. [5], then the result 
for 0(2) should have a higher systematic uncertainty. 
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FIG. 8. Comparison of NLO large N result with simulation results for 0(1), 0(2) and 0(4). The solid 
line is the NLO large N result for A{<j)'^)c/Nu, and the dashed line is the N=oo value. 
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0.008149(54) 
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8 
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2.5 
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144. 


12 
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32 


6. 
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3.5 
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288. 


24 
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-0.0000604(44) 


4.3 
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32 


12. 


0.00272961(71) 


-0.0002622(36) 
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48 


12. 
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-0.0004271(50) 


11.3 
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768. 


128 
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0.00213384(75) 


-0.0004142(63) 


1.6* 
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768. 
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32 
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3.7* 
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1536. 
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12. 
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7.5* 


1411 



TABLE III. Collection of 0(1) simulation data 



APPENDIX A: TABULATED DATA 



Tables III and IV are a collection of all the data reported in this paper. The standard 
integrated decorrelation time r for a single operator O is defined 



T 



1 ~ C{n 



E (Al) 



2 



where 



N-n 



i=l 




Cin) - E A.A^n - li^ E ^0 ' (A2) 



is the auto-correlation function associated with the operator A. In practice the sum in (Al) 
is cut-off when C(n)/C(0) first drops below 0.05 because of the statistical fluctuation in 
C(n). The nominal decorrelation time hsted in Tables III and IV is the largest value of the 
various operators required in the computations of the Binder cumulant and A(0^)c/w by 
using the canonical reweighting method ( [5]). 
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TABLE IV. Collection of 0(4) simulation data 
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